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Gamma rays induce a p53-independent mitochondrial 
biogenesis that is counter-regulated by HIF1a 

A Bartoletti-Stella 1 , E Mariani 1 , 1 Kurelac 1 , A Maresca 2 ' 3 , MF Caratozzolo 4 , L lommarini 5 , V Carelli 2 ' 3 , 
LH Eusebi 6 , A Guido 7 , G Cenacchi 8 , L Fuccio 6 , M Rugolo 5 , A Tullo 4 , AM Porcelli*' 5 ' 9 and G Gasparre*' 1 

Mitochondrial biogenesis is an orchestrated process that presides to the regulation of the organelles homeostasis within a cell. 
We show that y-rays, at doses commonly used in the radiation therapy for cancer treatment, induce an increase in mitochondrial 
mass and function, in response to a genotoxic stress that pushes cells into senescence, in the presence of a functional p53. 
Although the main effector of the response to y-rays is the p53-p21 axis, we demonstrated that mitochondrial biogenesis is only 
indirectly regulated by p53, whose activation triggers a murine double minute 2 (MDM2)-mediated hypoxia-inducible factor 1a 
(HIF1a) degradation, leading to the release of peroxisome-proliferator activated receptor gamma co-activator 1/f inhibition by 
HIF1a, thus promoting mitochondrial biogenesis. Mimicking hypoxia by HIF1a stabilization, in fact, blunts the mitochondrial 
response to y-rays as well as the induction of p21 -mediated cell senescence, indicating prevalence of the hypoxic over the 
genotoxic response. Finally, we also show in vivo that post-radiotherapy mitochondrial DNA copy number increase well 
correlates with lack of HIF1 a increase in the tissue, concluding this may be a useful molecular tool to infer the trigger of a hypoxic 
response during radiotherapy, which may lead to failure of activation of cell senescence. 
Cell Death and Disease (2013) 4, e663; doi:1 0.1 038/cddis.201 3.1 87; published online 13 June 2013 
Subject Category: Cancer 



Mitochondrial mass and function are controlled through 
biogenesis of mitochondria, which is pivotal to maintain 
energetic homeostasis. Peroxisome-proliferator activated 
receptor gamma co-activators (PGC) preside to such process 
via a cascade of nuclear proteins that ultimately regulate the 
balance between the organelles fusion and fission, the 
replication and transcription of the mitochondrial genome 
(mtDNA), the assembly and function of the respiratory chain 
components and the homeostasis of all metabolic enzymes 
that reside within mitochondria. 1-3 Such regulation is carried 
out in a two-way manner and is an integrated pathway that 
requires a cross-talk between nucleus and mitochondrion, 4,5 
since the latter is a semi-autonomous organelle whose 
replication appears to be independent from the cell cycle, 6-8 
and may follow its own checkpoints. 9 Mitochondria are 



considered as the principal oxygen sensor of the cell, 10 as 
one of their main roles is adenosine tri-phosphate (ATP) 
production via oxidative phosphorylation, which requires 0 2 
as the final acceptor of the electron transfer. Another pivotal 
oxygen sensor is the master regulator of the hypoxic 
response, namely hypoxia-inducible factor 1a (HIF1a), a 
ubiquitous transcription factor that is kept, under normoxia, 
at often undetectable protein level by continuous degradation 
via the 26S proteasome. 11 Cross-talk between the two 
sensors is necessary to shut off oxidative phosphorylation 
and shift metabolism towards oxygen-independent glycoly- 
sis 12 ' 13 in shortage of oxygen. Activation of hypoxia-induced 
pathways is largely investigated in cancer, where it is thought 
to strongly contribute to a failure in responding to y-radiation 
therapy. 14 Gamma rays have been long exploited as an 
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anticancer strategy due to their ability to cause double-strand 
breaks within DNA and induce a genotoxic response that locks 
the cell in a non-replicative state to allow repair. 1516 The 
canonical DNA damage response (DDR) is orchestrated by 
several protein kinases transducing stress signals, whose 
activation propagates the damage signal to effector molecules, 
the most important of which being the tumour suppressor 
p53. 17 ' 18 N-terminal phosphorylated p53 traslocates to the 
nucleus, escaping a constitutive degradation mediated by the 
murine double minute human orthologue M DM2. 1920 p53 
favours replication arrest via, among others, the transcription of 
cell-cycle regulator cyclin-dependent kinase inhibitor 1A 
(CDKN1A, alias p21). Such p53-p21 axis is one of the major 
pathways presiding to the trigger of cell senescence, the 
induction of which represents at the same time the goal and the 
drawback of anticancer radiotherapy. 21 In fact, senescent cells 
are metabolically active yet quiescent, acquire widespread 
morphological changes and become resistant to apoptosis. 
Hence, cancer cells in a senescent state might remain as 
'dormant' malignant cells, and therefore represent a dangerous 
potential for tumour relapses, since through disruption of p53 
and p21 function they may escape senescence. 22 ' 23 None- 
theless, radiotherapy is one of the most effective non-surgical 
treatments of cancer patients with better functional preserva- 
tion and less systemic influences. 24 The understanding of the 
complex interplay between genotoxic and hypoxic responses 
remains a critical goal to predict therapy outcome. In fact, 
prevalence of one over the other is thought to determine the 
fate of cancer cells in response to therapy, while the cross- 
regulation of the mechanisms underlying activation, stabiliza- 
tion, and degradation of p53 and HIF1a, and their interaction 
with MDM2, remain largely to be clarified 25 At the same time, 
identification of molecular markers that may function as 
predictors of y-rays therapy outcome may be useful prognostic 
tools. In this context, y-radiation has been shown to induce an 
increase in mtDNA copy number in the brain and spleen tissue 
of mice. 26 Similarly, X-rays were shown to induce increase in 
mitochondrial mass and oxidative function under the cell-cycle 
control, 27 whereas oncocytic relapses, tumours with an 
aberrant mitochondrial biogenesis 28 were observed to occur 
post-radiotherapy of colorectal cancers. 29 ' 30 Overall, these 
findings point at mitochondrial homeostasis as one of the 
biological processes most affected by ionizing radiations. 

We here investigate the complex mechanistic relationships 
between p53-genotoxic and HIF1a-hypoxic responses in the 
regulation of mitochondrial biogenesis and cell senescence 
following y-rays treatment. We demonstrate that y-rays strongly 
induce mitochondrial biogenesis, and that the entity of such 
process is dependent on a p53-HIF1 a regulated axis that relies 
on MDM2-mediated HIF1a degradation. Last, we attempt to 
translate our findings in vivo and set the basis for the use of 
mitochondrial parameters to infer efficacy of radiotherapy. 

Results 

y-Radiation induces a net mitochondrial biogenesis. To 

test the effect of y-rays on mitochondrial biogenesis, we 
evaluated the mtDNA abundance in irradiated cells. We first 
analysed TP53 wild-type cells, namely RPE1 and HCT116. 
As expected, p53 was shown to significantly increase upon 



radiation treatment, indicating activation of a genotoxic 
response, which was further confirmed upon increased 
expression of the p53-dependent p21 (Figure 1a; 
Supplementary Figure 1A). Both RPE1 and HCT116 cells 
showed mtDNA copy number increase both at the last day of 
irradiation (IR) and after 120h recovery (REC) in culture 
(Figure 1b), suggesting that the copy number increase is an 
ongoing process proceeding beyond stimulus withdrawal. In 
agreement with the coupling between mtDNA replication and 
transcription 31 the latter was also observed to increase 
following IR at 120h from stimulus withdrawal (Figure 1b). 
Mitochondrial proteins were also shown to have an increased 
expression upon IR, consistently with a complete mitochon- 
drial biogenesis activation. In fact, proteins such as the fusion 
regulator optic atrophy 1 (OPA1) and the mitochondrial 
transcription factor A (TFAM), involved in mtDNA mainte- 
nance and transcription, showed increased expression upon 
IR and at recovery. The same trend was observed for 
voltage-dependent anion channel 1 (VDAC1), an indicator of 
mitochondrial mass. Similarly, expression of mitochondrial 
respiratory complexes subunits (complex V ATP5/?, complex 

II (CM) succinate dehydrogenase subunit A/B (SDHA and 
SDHB), complex I (CI) NDUFA9 and NDUFS3, and complex 

III Core II) was shown to increase in response to IR 
(Figure 1c; Supplementary Figure 1B). Gene expression 
evaluation of peroxisome-proliferator activated receptor 
gamma co-activator ^p (PGC-1/?), one of the master 
regulators of mitochondrial biogenesis, indicated a radia- 
tion-dependent increase in both RPE1 and HCT116 up to 
two- and three-fold, respectively (Figure 1d). Expression of 
PGC-1a, the other biogenesis co-activator, was undetectable 
in both RPE1 and HCT116, consistently with the knowledge 
that it is expressed at high levels in tissues where oxidative 
metabolism is particularly active, such as brown adipose 
tissue or skeletal muscle. 32 Concordantly, downstream PGC- 

target genes were also shown to be upregulated upon IR 
(Supplementary Figure 2). These data suggested an overall 
increase in mitochondrial mass (Figure 1e), which was 
confirmed by fluorescence-activated cell sorting (FACS) 
analysis (Supplementary Figure 3A and B), likely to be 
driven by a PGC-1/? increase. To understand whether this 
phenomenon was due to a compensatory effect following 
damage, we evaluated the bioenergetics competence of both 
cell lines. Following IR, both CI- and Cll-driven ATP 
synthesis were shown to be significantly increased over 
two-fold in RPE1 and likewise in HCT116 (Figure 1f). The 
same result was obtained upon evaluation of the activity of 
citrate synthase (CS), an enzymatic marker of mitochondrial 
mass (Figure 1f). However, when ATP synthesis and CS 
activity were normalized over protein content, no changes 
were observed (Figure 1g), suggesting that the induced 
mitochondrial biogenesis is more likely due to a coordinated 
cell and mitochondrial mass increase than to a compensatory 
mechanism for some radiation-induced mitochondrial 
damage, and that cells remained metabolically active. 

y-Rays-induced mitochondrial biogenesis is indepen- 
dent on p53. Since PGC-1 p appeared to be the main driver 
of the y-rays-induced mitochondrial biogenesis, and following 
y-IR p53 is activated, we sought to determine a potential role 
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Figure 1 Gamma rays induce mitochondrial biogenesis in RPE1 and HCT1 1 6. (a) Western blot analysis of p53 and p21 at the last day of irradiation (IR) and after 120 h of 
recovery (REC). Gamma tubulin was used as a loading control. One representative experiment of three is shown, (b) Relative mtDNA copy number and mtRNA transcription. 
Data are mean ± S.D. (n = 3, ANOVA, *P< 0.05 versusJO). (c) Western blot analysis of mitochondrial proteins at TO, IR, and REC. (d) Relative mRNA quantification of PGC- 
1 p upon IR. Data are mean ± S.D. (n = 3, f-test, *P< 0.05 versusJO). (e) Cells stained with MitoTracker Red (MTR) to evaluate mitochondrial mass and network morphology. 
Calcein-AM staining was used to highlight cells morphology. One representative experiment of three is shown. Magnification x 63/1.4. Length bar corresponds to 25 ^m. 
(f) The rate of mitochondrial ATP synthesis driven by complex I (CI) and complex II (CI I) substrates, and citrate synthase (CS) activity after irradiation. Data are mean ± S.D. 
(n=3, f-test, *P<0.05 versusJO). (g) ATP synthesis driven by CI and Cll and CS activity normalized on protein content. Data are mean ± S.D. (n = 3) 



for p53 in the phenomena we observed. p53 has been 
recently shown to display a repressor activity at the 
promoters of murine pgc-1u and -ip genes, albeit in the 
context of telomere dysfunction, so that overexpressed p53 
binds its specific responsive elements (REs) to reduce 
expression of both regulators. 33 Such an analysis has not 
been carried out on the orthologue human genes. We looked 
for specific p53 REs in the promoter of human PGC-1 ft, and 
detected two close decamers starting at position -2732; 
-2706 (GRCh37-EntrezGene ID: 133522; Figures 2a and 
b). We hence took advantage of the syngenic HCT116 cell 
line knock-out (KO) for TP53, 34 namely HCT1 16 TP53K °. 
Complementation with wild-type p53 (Figure 2c) did not 
prevent PGC-1 expression increase following IR, similarly 
to what occurred in mock HCT116 KO cells (Figure 2d). 
Consistently, p53 was shown not to bind to its PGC-1 /? weak 
REs (Figure 2e). p53 activation was confirmed by increase in 
p21 expression upon y-IR (Supplementary Figure 4) and 
through the binding of p53 to its REs in the p21 promoter 
(Figure 2e). 

To rule out a potential role for p53 activation in triggering the 
mitochondrial effects induced by y-rays, we evaluated 



mitochondrial biogenesis parameters in HCT116 T 



and 



in an osteosarcoma-derived cell line, HPS11, whose TP53 
harbours a hemizygous mutation within the DNA binding 
domain, that is, the c.13055G>C, p.156R>P (Figures 3a 
and b). Complementation of HCT1 16 TP53KO cells with a 
p.156R>P-mutated p53 confirmed the lack of p21 transacti- 
vating ability of such mutant compared with HCT116 TP53KO 
complemented with wild-type p53 (Figure 3c), as previously 
observed. 35 

Both HPS11 and HCT116 TP53KO were more sensitive to 
radiation treatment compared with RPE1 and HCT116, not 
surviving beyond the fourth dose of radiation (Supplementary 
Figure 5). Nonetheless, both HCT116 TP53KO and HPS11 
displayed a 2.5-fold y-rays-induced mtDNA copy number 
increase, similar to that occurring in RPE1 and HCT116, as 
well as an increase in mtDNA transcription (Figure 3d). 
Moreover, levels of PGC-1 ft transcript and expression of 
mitochondrial proteins, both structural and belonging to 
oxidative phosphorylation, were shown to markedly increase 
upon IR (Figures 3e and f; Supplementary Figure 6), 
concordantly with cellular and mitochondrial mass increase 
(Figures 3g and h; Supplementary Figure 3C). 
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Figure 2 p53 is not a PGC-1 p regulator, (a) Schematic representation of the genomic location of the p53-responsive elements (REs) within the promoter region of PGC- 
1/?. The genomic coordinates are relative to the human genome build hg19. (b) Promoter region of human PGC-1/?, stretching from -2845 to -2542 with respect to 
transcription starting site (TSS). The putative p53-RE is underlined, (c) Western blot analysis of p53 after HCT1 16 TP53KO complementation with wild-type p53. Gamma tubulin 
was used as a loading control. One representative experiment of three is shown, (d) Expression of PGC-1/? both in HCT116 TP53KO and in HCT1 1 6 p53 upon IR. Data are 
mean ± S.D. (n = 3, f-test, *P<0.05 versus TO), (e) p53 binding and histone H4 acetylation at the promoter of PGC-1 y5 in RPE1 cells upon IR. Data are mean ± S.D. (n = 2, 
f-test, *P<0.05). The occupancy of p53 and H4 acetylation at the promoter of p21 in the same cells is used as positive control 



Altogether, these findings point to a lack of involvement of 
p53 in orchestrating the cell response to y-rays in terms of 
activation of mitochondrial biogenesis. 

HIF1a stabilization blunts y-rays-induced mitochondrial 
biogenesis activation. Next, we sought to determine 
whether p53 activation may influence effectors that may 
instead inhibit mitochondrial biogenesis, such as HIF1a. In 
fact, it is well known that genotoxic and hypoxic stress 
responses are interwoven and tightly regulated by proteins 
such as MDM2. p53 and HIF1a may be sequestered or 
activated, or even led to degradation by MDM2, according to 
the prevalence of one or the other response. 25 The relation- 
ship between HIF1a and mitochondrial biogenesis has been 
previously detailed, especially with respect to its negative 
PGC-1/? regulation via the c-MYC oncogene. 36 Following IR, 
HCT116 cells were shown to have a time-dependent HIF1a 
protein decrease, which was most profound at recovery 
(Figure 4a; Supplementary Figure 7A), and therefore 
antithetically followed p53 protein expression (Figure 1a). 
Such decrease was a post-transcriptional effect since HIF1a 
mRNA was shown to increase at recovery (Figure 4b). 
Interestingly, an MDM2-dependent HIF1a degradation has 
been previously reported. 37 To validate such a mechanism in 
the cells analysed, based on the observation that MDM2 also 
increased following IR (Figure 4c; Supplementary Figure 7B), 
we overexpressed MDM2 and followed HIF1a degradation 



after blocking protein translation with cycloheximide (chx). 
HIF1a was shown to be rapidly degraded only when MDM2 
was in excess within the cells, even without a radiation 
stimulus, indicating that a basal MDM2-dependent HIF1a 
degradation occurs in HCT116 (Figure 4d; Supplementary 
Figure 8A). Such degradation did not occur when MDM2 was 
overexpressed in HCT1 16 that also overexpressed an HIF1a 
mutated on the three residues responsible for hydroxylation- 
dependent degradation, 38 suggesting these residues to be 
essential (Figure 4e; Supplementary Figure 8B). The 
necessity for hydroxylases to mediate MDM2-dependent 
degradation was further confirmed when the potent HIF1a 
stabilizer dimethyloxallylglycine (DMOG), a specific inhibitor 
of prolyl-hydroxylases (PHDs) activity, was used. In such 
conditions, HIF1a protein did not decrease upon chx 
treatment (Figure 4f; Supplementary Figure 8C). We hence 
attempted to force a hypoxic response concomitantly to the 
genotoxic one, with the aim to reveal whether mitochondrial 
biogenesis would be activated according to the prevalence of 
one or the other. DMOG-treated HCT1 1 6 cells were shown to 
retain stabilized HIF1a for the entire duration of the radiation 
treatment along with increased expression of canonical 
HIF1a-responsive genes lactate dehydrogenase A (LDHA) 
and BCL2/adenovirus E1B 19KDa interacting protein 3-like 
(BNIP3L) (Supplementary Figure 9), albeit not displaying any 
increase in mtDNA copy number (Figure 4g). The same 
finding was obtained when cells were treated with HIF1a 
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Figure 3 p53 does not directly regulate mitochondrial biogenesis induced by y-rays. (a) ElectropherogramsoftheHPSH cells inactivating 7P53 mutation (p. 156R>P) within 
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f-test, *P< 0.05 versus TO), (e) Relative mRNA quantification of PGC-1 p in HPS1 1 and HCT1 1 6 TP53KO cells after irradiation. Data are mean ± S.D. (n = 3, f-test, *P< 0.05 versus 
TO), (f) Western blot analysis of mitochondrial proteins upon IR. (g and h) MitoTracker Red (MTR) and Calcein-AM staining of HPS1 1 (g) and HCT1 16 TP53KO (h) cells to evaluate 
mitochondrial mass and network morphology prior (TO) and after irradiation (REC). Magnification x 63/1 .4. Length bar corresponds to 25 



stabilizer desferrioxamine (DFO) (Supplementary Figure 10). 
Concordantly, neither an increase in PGC-1 ft expression nor 
of its target genes (Supplementary Figure 11) after IR was 
observed in DMOG-treated cells (Figure 4h). Neither mtDNA 
maintenance proteins nor respiratory complexes subunits 
were shown to increase after IR in DMOG-treated cells 
displaying a chronically stabilized HIF1a (Figure 4i; 
Supplementary Figure 12). Overall, no increase in mitochon- 
drial mass was observed, despite a net cell mass increase 
(Figure 4j). Analogous findings were obtained when RPE1 
cells were treated with DMOG and irradiated (Supplementary 
Figure 13). Overall, these data indicate that the hypoxic 
response prevails over the genotoxic one in terms of 
mitochondrial biogenesis activation that may be induced by 
IR via an indirect mechanism, namely a re-positioning of 
MDM2 from p53 to HIF1a. This would foster hydroxylated 
HIF1a degradation upon p53 activation, preventing HIF1a to 
inhibit mitochondrial biogenesis. 

HIF1a prevents radiation-induced cell senescence. We 

next attempted to investigate whether, along with mitochon- 
drial biogenesis, the hypoxic response would prevail in terms 
of repressing y-rays-induced cell senescence, since this 
would allow to link the two biological processes. 



Mitochondrial biogenesis parameters would then be useful 
to monitor the efficacy of I R in inducing senescence, which is 
the main scope of radiotherapy. Interestingly, the stabiliza- 
tion of p53, its main targets p21 and MDM2 did not occur 
after IR in DMOG-treated RPE1 and HCT116 cells, nor did 
p16 increase, another senescence effector (Figure 5a; 
Supplementary Figure 14). The genotoxic stimulus was not 
able to activate at all p53 in the presence of stabilized HIF1a, 
as no binding of p53 could be detected on the p21 promoter 
in DMOG-treated cells, and consequently no H4 histone 
acetylation occurred (Figure 5b). Accordingly, the expression 
of activated p53 gene targets BAX and GADD45A was not 
shown to increase, indicating the lack of p53 transcriptional 
activity (Figures 5c and d). 

An altered morphology biased towards a cell enlargement, 
along with the proliferative arrest, but not death, of the cell 
population following the optimized Gray doses, are clear 
hallmarks of cell senescence. 39 ' 40 HIF1a stabilization by 
DMOG was able to prevent cells to enter senescence, as 
DMOG-treated cells did not undergo morphological changes 
and did not show a different beta-galactosidase (/?-Gal) 
staining compared with controls after IR (Figures 5e and f). 
The same effect could be observed when DFO was used 
(Supplementary Figure 15). Taken together, these data 
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Mest, *P<0.05 versus TO), (c) Western blot analysis of MDM2 in HCT1 16 cells after irradiation; y-tubulin was used as a loading control. One representative experiment of 
three is shown, (d) Western blot analysis of HIF1a in HCT116 cells overexpressing HIF1a (HCT116 HIF1 a) upon overexpression of MDM2 ( + ) compared with clones not 
overexpressing MDM2 ( - ). Degradation of HIF1 a was evaluated by treating cells with 50 /^g/ml of cycloheximide (chx) for 0, 1 0, and 20 min; y-tubulin was used as a loading 
control, (e) Degradation evaluation of HIF1a mutated on the three residues subjected to hydroxylation (HIF1a), in the presence and absence of overexpressed MDM2. 
Degradation of HIF1 a was evaluated by treating cells with 50 ^g/ml of chx for 0, 10, and 20 min; y-tubulin was used as a loading control, (f) Degradation evaluation of HI F1 a in 
the presence and absence of overexpressed MDM2, upon forced pharmacological HIF1 a stabilization by 1 DMOG. Degradation of HIF1a was evaluated by treating cells 
with 50 fig/ml of chx for 0, 10, and 20 min; y-tubulin was used as a loading control, (g) Relative mtDNA copy number and (h) PGC-1 p expression after irradiation in HCT1 1 6 
cells with stabilized HIF1a after 1 DMOG treatment. Data are mean ± S.D. (n=3, Mest, *P<0.05 versus TO), (i) Western blot analysis of mitochondrial proteins in 
HCT116 cells treated with DMOG and irradiated, (j) Staining with MitoTracker Red (MTR) and Calcein-AM to evaluate mitochondrial mass and network morphology. 
Magnification x 63/1 .4. Length bar corresponds to 25 ^m 



indicate that, along with mitochondrial biogenesis, HIF1a 
stabilization is able to blunt the cell senescence response, 
suggesting the two processes to be concomitant. 

Last, we tested the hypothesis that y-rays may trigger 
mitochondrial biogenesis in wVoand that an easy-to-measure 
parameter, such as mtDNA copy number, might be predictive 
of HIF1a stabilization. We took advantage of a collection of 
healthy tissue rectal biopsies taken from subjects who had 
undergone radiation treatment for prostate cancer, since y- 
rays also hit the rectum lining. We evaluated mtDNA copy 
number in rectal epithelium before IR compared with biopsies 
taken 2 months after the end of therapy. Surprisingly, mtDNA 
copy number was shown to be significantly higher in 16/32 
(50%) subjects (Figure 6a). Considering all data together, the 
increase was significant over all subjects analysed (pre-IR 
versus post-IR, P<0.05), suggesting that the radiation trigger 
had lasted long after the withdrawal of the stimulus 
(Figure 6a). A significant inverse correlation was found 
between decrease in HIF1a protein and mtDNA copy number 
increase (Figures 6b and c; Fisher's exact test P= 0.039), 
strengthening the negative relationship between HIF1a 
stabilization and y-ray-induced mitochondrial biogenesis. 



Discussion 

In this paper we show that y-rays induce a p53-independent 
net mitochondrial biogenesis blunted by HIF1a stabilization. 
An increase in mitochondrial mass and function had been 
previously observed in lung carcinoma cells after exposure to 
X-rays. 27 Our results extended such findings to y-rays and 
thoroughly characterized the mitochondrial changes that 
occur when a genotoxic stimulus is provided. In fact y-rays 
elicit a DDR that, in presence of a functional p53, compels 
cells to undergo a temporary replicative arrest until repair. 15,40 
One of the critical mediators of such arrest is p21, a 
transcriptional target of p53, which binds to and inhibits 
Cdk2/cyclin E complexes, thereby delaying the G1/S phase. 41 
Unsurprisingly, p21 also has a pivotal role in triggering and 
maintaining cell senescence, the permanent replicative block 
cells undergo upon entering a quiescent state while remaining 
metabolically active, which is the main achievement of ir- 
radiation when treating cancer. 21 Despite being regulated 
mainly through activation of p53, p21 expression appears to 
be governed also by other factors, such as E2A and N-myc 
downstream regulated gene 1 (NDRG1), 42 ' 43 which may 
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Figure 5 HIF1 a stabilization blunts cells genotoxic response and cell senescence, (a) Western blot analysis of p53, p21 , MDM2, and p16 in RPE1 and HCT1 16 cells in the 
absence and presence of 1 [M DMOG and irradiated. Gamma tubulin was used as a loading control. One representative experiment of three is shown, (b) p53 occupancy and 
histone H4 acetylation at the promoter of p21 in irradiated RPE1 cells, in the presence or absence of stabilized HIF1a by 1 /uM DMOG. Data are mean ± S.D. (n=2, f-test, 
*P<0.05 versus TO), (c and d) Relative expression of p53-responsive genes BAXand GADD45A in RPE1 cells, in the presence or absence of stabilized HIF1a by 1 ^M 
DMOG, and irradiated. Data are mean ± S.D. (n = 3, f-test, *P< 0.05 versusJO). (e and f) Senescence-associated /?-galactosidase staining in RPE1 and HCT1 16 cells in the 
absence and presence of 1 jM DMOG and irradiated. One representative experiment of three is shown. Magnification x 10 



explain p21 expression changes in the absence of p53 
activation (Figure 5a). 

Replication and function of mitochondria are semi-autono- 
mous processes. Nuclear-encoded proteins are needed for 
the organelles biogenesis, such as the DNA polymerase 
gamma or OPA1, yet checkpoints for mtDNA replication, for 
instance, do not seem to exist or are hypothesized to be 
different from those of nuclear DNA. 9 The increase in 
mitochondrial mass and function we observed in irradiated 
cells paralleled the cell mass increase that is a hallmark of cell 
senescence. 39 Such results find explanation in the semi- 
autonomous replicative mode of mitochondria, which would 
therefore continue to proliferate despite lack of cell division, in 
the senescent yet metabolically active G1 -arrested cells. The 
mitochondrial biogenesis increase observed after y-rays 
treatment, therefore, may not be a compensatory effect since 
mitochondrial function increases in irradiated cells. Albeit the 
increase in mitochondrial biogenesis appeared independent 
of p53 activation or even presence, unlike in mice, 33 the 



influence of p53 on mitochondrial biogenesis may be indirect, 
via the known MDM2-mediated HIF1a degradation, 37 a 
hypothesis supported by the increase in MDM2 observed in 
HCT116 TP53KO and in HPS11 cells (Supplementary 
Figure 16). Post-radiation MDM2 increase may hence also 
occur independently of p53 stabilization, either through other 
members of the p53 family, such as p73, 44 or through an 
NFjcB-mediated induction, 45 and contribute to HIF1a degra- 
dation and mitochondrial biogenesis increase we observed in 
p53 KO cells. Both HIF1a and the mitochondrial respiratory 
chain are considered as the two fundamental oxygen sensors 
of the cell. As 0 2 tension decreases, HIF1a stops becoming 
hydroxylated and is hence stabilized and able to transcribe a 
large set of genes involved in the regulation of glycolysis and 
neoangiogenesis 46 47 Since 0 2 is also the major substrate for 
the electrons of the respiratory chain, the latter is slowed down 
during hypoxia, and cells must rely on non-oxidative metabo- 
lism for ATP production. A basal HIF1 a degradation by MDM2 
occurs in cells, 37 with mechanisms that are likely to be 
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Figure 6 Counter-regulation of HIF1 a expression and mtDNA copy number after radiation treatment in vivo, (a) Absolute mtDNA copy number in rectal biopsies of patients 
subjected to radiotherapy for prostate carcinoma, prior to (pre-IR) and post-irradiation (post-IR). Data are mean ± S.E.M. (*P<0.05; paired f-test), Wilcoxon Rank tests were 
used for median comparison of all pre- and post-irradiation biopsies (*P<0.05 versus Pre-IR). (b) Western blot analysis of the expression of HIF1a in representative cases in 
which the mtDNA copy number was found not to change or decreased after radiation treatment, (c) Western blot analysis of the expression of HIF1 a in representative cases in 
which mtDNA copy number increased after radiation treatment. Beta-actin (ACT) was used as a loading control. One representative blot of three is shown. Fisher's exact test 
was conducted to evaluate significance of the correlation between HIF1a protein level increase (calculated based on densitometric analysis as a ratio between post-radiation 
and pre-radiation protein level greater than 1, in all samples) upon irradiation and lack of mtDNA copy number increase 



analogous to the 0 2 -dependent hydroxylation by PHDs, since 
our mutant HIF1 a stops responding to MDM2 degradation. It is 
therefore plausible to envision that when p53 is activated upon 
a genotoxic stimulus, such as y-rays, MDM2 may not preside 
to its degradation and be re-directed towards HIF1a, whose 
degradation would hence be fostered (Figure 7). The relation- 
ship between HIF1a and mitochondrial biogenesis has been 
previously investigated, given the fact that the two processes 
are obviously antithetical. It has been shown that HIF1a is an 
indirect repressor of PGC-1/? through a c-MYC-mediated 
axis, 36 which overall leads to a downregulation of mitochon- 
drial biogenesis during hypoxia. Our findings therefore point to 
an indirect control of mitochondrial biogenesis by activated 
p53, whose role would be to allow MDM2 to degrade HIF1a, 
with a subsequent release of the basal inhibition by the latter 
over PGC-1/?. Our data also shed light on the controversial 



role of MDM2 in orchestrating the balance between genotoxic 
and hypoxic response. 25 ' 48 Although it has been suggested 
that MDM2 may either degrade or stabilize HIF1a, our results 
support the former hypothesis, at least upon y-rays treatment. 

A rather important finding is that when hypoxia is induced in 
parallel to a genotoxic stimulus, the former response prevails 
over the latter. In fact, cells treated with DMOG are prevented 
to undergo senescence when irradiated, as indicated by the 
unequivocal /?-Gal marker and by the lack of p21 increase. 
Interestingly, p53 stabilization also appears to be prevented, 
likely since non-hydroxylated HIF1a may not be recognized 
for degradation by MDM2, which would then re-direct over to 
p53. Further investigation is warranted to understand in detail 
why the p53-p21 axis appears to be repressed upon HIF1a 
stabilization. In terms of mitochondria, which are the focus of 
our paper, it is most interesting that the y-rays-induced 
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Figure 7 Schematic representation of the cross-regulation of genotoxic and hypoxic response upon gamma irradiation. Under basal conditions, MDM2 binds and 
degrades p53 (1). In normoxia, upon gamma radiation stimulus (2), p53 activation promotes the transcription of p21 and MDM2an6 allows MDM2 to detach from p53 and 
redirect over to HIF1 a (3). In the absence of p53, the accumulation of MDM2 following gamma radiation may also be due to other factors such as p73 and NF-kB (4), however 
leading to HIF1a degradation (5). The HIF1a-mediated repression of PGC-1/? is therefore inhibited (6), which enables activation of mitochondrial biogenesis (7). Cell 
senescence triggered by the activated p53/p21 axis (8), may instead not occur under hypoxia, when HIF1 a is stabilized and MDM2 is redirected to p53. Functionally and non- 
functional^ active factors are marked in red and grey, respectively 



biogenesis is also strongly blunted, in parallel to senescence, 
during hypoxia. We have managed to show this to occur 
in vivo as well, 2 months after the withdrawal of the y-rays 
stimulus. Tissues hit by y-rays retained a 'scar' of activated 
mitochondrial biogenesis only when local hypoxia had not set 
in, as shown by the decreased HIF1 a expression. Whether the 
lack of a net mitochondrial biogenesis increase is due to a 
continuous cell replication that does not allow organelles to 
accumulate, or to a yet unknown and more direct molecular 
mechanism, biogenesis parallels the activation of cell 
senescence and contributes to the signalling of an ongoing 
genotoxic response. In this context, triggering of a hypoxic 
response may compromise efficacy of y-rays in inducing 
senescence, as we have here shown. Evaluation of mito- 
chondrial biogenesis parameters in vivo, such as mtDNA copy 
number, may therefore reveal to be an efficient tool to predict 
the trigger of such a response. 

Materials and Methods 

Cell cultures. Human colorectal carcinoma cell line HCT116, their genetically 
modified derivates lacking TP53, 34 HCT1 1 6 TP53KO , and the osteosarcoma cell line 
HPS11 were grown in Dulbecco's Modified Eagle's Medium (DMEM) supple- 
mented with 10% fetal bovine serum (FBS), 100 units/ml penicillin, 100 ^g/ml 
G-streptomycin and 2 mmol/l L-glutamine. Non-tumour immortalized RPE1 (human 
retinal epithelial) cells were cultured in DF12 medium (DMEM/nutrient mixture F-12 
ham; Sigma-Aldrich, Milan, Italy) supplemented with 10% FBS and 100 units/ml 



penicillin, 100/zg/ml G-streptomycin. Cell lines were maintained in a humidified 
incubator at 37°C with 5% C0 2 . 

y-Rays treatment and HIFta compound-mediated stabilization. 

Treatment was optimized to (i) use y-rays doses close to those currently used in 
radiotherapy, (ii) maintain cells viable (Supplementary Figure 2), and (iii) eliminate 
dead cells from the population. In all, 200000 cells were seeded per 25 cm 2 flask 
(TO). The following day flasks were submitted to a dose of 4 Gray (Gy), specifically 
2Gy/min with the biological irradiator IBL437C (89-294). The treatment was 
repeated after exactly 24 h for 4 days (IR). Alternatively, cells were allowed to 
recovery for 120 h before harvesting (REC). HIF1a was stabilized by adding to the 
culture medium the inhibitor of proline hydroxylase DMOG (1 /M) or the iron 
chelator DFO (250 nM) (Sigma-Aldrich) for the entire duration of the radiation 
treatment. 

Western blotting. Total lysates (80 ^g) were separated by 10-12% sodium 
dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE), and nitrocellu- 
lose membranes were incubated with primary antibodies against y-tubulin (Sigma- 
Aldrich), HIF1a (Bethyl Laboratories, Montgomery, TX, USA), p53, p21, p16 
(Santa Cruz Biotechnologies, Dallas, TX, USA), OPA1 (Bioscience, Oxford, UK), 
TFAM (AbCam, Cambridge, UK), SDHA, SDHB, NDUFA9, NDUFS3 (Mito- 
Sciences, Eugene, OR, USA), ATP5/? (BioVision, Milpitas, CA, USA), VDAC 
(Millipore, Billerica, MA, USA), Core II (Life Technologies, Monza, Italy), and 
MDM2 (Calbiochem, Billerica, MA, USA). Densitometric analyses were performed 
with Image J software. 

Cellular and mitochondrial morphology. Cellular and mitochondrial 
morphology was assessed after cell staining with 100nM calcein-acetoxymethyl 
ester (calcein-AM, Life Technologies) and 10nM MitoTracker Red (MTR, Life 
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Technologies) for 30 min at 37°C. Images were acquired with a digital imaging 
system using an inverted epifluorescence microscope with x 63/1.4 oil objective 
(Diaphot, Nikon, Japan) and a back-illuminated Photometries Cascade CCD 
camera system (Roper Scientific, Tucson, AZ, USA). Images were captured and 
analysed using the Metamorph software (Universal Imaging Corp., Downingtown, 
PA, USA). 

ATP synthesis evaluation. The rate of mitochondrial ATP synthesis driven 
by CI and CI I was performed in digitonin-permeabilized cells and normalized on 
protein content, as previously described. 49 Briefly, aliquots of cells (0.1-0.2 mg 
protein) were incubated with 5 mM malate plus 5 mM pyruvate (CI substrates) or 
with 10 mM succinate (Cll substrate) plus 2^g/ml rotenone. The reaction was 
started by addition of 0.2 mM ADP in the presence of luciferine/luciferase, and 
chemiluminescence was evaluated as a function of time with a luminometer. After 
addition of 10mM oligomycin, the chemiluminescence signal was calibrated with 
an internal ATP standard. CS activity was determined on an aliquot of 
permeabilized cells, as previously described. 50 

mtDNA copy number. mtDNA content was evaluated by a Real Time-PCR 
multiplex assay based on hydrolysis probe chemistry. Primers, probe, and 
conditions were previously published. 51 About 2 ng of DNA extracted from cells 
and tissues was analysed using the LightCycler 480 Probes Master Mix (Roche 
Diagnostics, Milan, Italy) and the LightCycler 480 (Roche Diagnostics) instrument. 

Gene expression analyses via real-time PCR. RNA was extracted 
from 10 6 cells with TRIZOL reagent (Life Technologies) following manufacturer's 
instructions. Quantitative real-time PCR (qRT-PCR) for MT-ND5, p21, GADD45A, 
BAX, PGC-10, HIF1a, POLG, MFN1, MFN2, COX5B, UQCRQ, NDUFB5, and 
TUBB1 was performed on complementary DNA (cDNA) (1 : 50) prepared by 
retrotranscription of 1 ^g of total RNA with High Capacity cDNA Reverse 
Transcription Kit (Life Technologies). The qRT-PCR was performed with GoTaq 
qPCR Master Mix (Promega, Madison, Wl, USA) and run in 7500 Fast Real-Time 
PCR System (Life Technologies), using following conditions: 95°C 5 min; 45 cycles 
of 95°C 15 s and 60°C 45 s. The calculations were performed following 2~ AACT 
method. Primer sequences were designed using Primer3 software. 52 Primer 
sequences and concentrations used in the real-time PCR are available on request. 

In silico promoter analyses. The promoter and 5'untranslated region 
(UTR) regulatory regions of the PGC-1/? gene were scanned to verify the 
presence of one or more p53 REs (p53REs) using the PatSearch algorithm 
implemented in the DNAfan tool. 53,54 The syntax pattern of the p53RE is made up 
of two tandem repeated decamers complying with a specific consensus 
corresponding to the 5'-PuPuPuC(A/T)(T/A)GPyPyPy-3' sequence, allowing at 
most three mismatches. The decamers can be spaced by 0-1 3 nt. The in silico 
analysis revealed the presence of one p53RE at position -2732; -2706 
(promoter region) upstream the transcription start site (TSS) of the PGC-1/? gene 
(GRCh37 - Entrez Gene ID: 133522). 

Chromatin immunoprecipitation. After treatments, proteins were cross- 
linked to DNA in living nuclei by adding formaldehyde directly to the cell culture 
medium to a final concentration of 1%. Cross-linking was allowed to proceed for 
10 min at 37°C. Cells were washed with phosphate-buffered saline, scraped off 
the plates and resuspended in 20 mM Tris-chloride pH 8.3, 3mM MgCI 2 , 20 mM 
KCI, 1 mM phenylmethylsulphonyl fluoride, 1 /uglml leupeptin, and 1 ^g/ml 
aprotinin. Nuclei were pelleted by microcentrifugation and lysed by incubation in 
nuclear lysis buffer (1% sodium dodecyl sulphate, 10mM EDTA, 50 mM Tris- 
chloride pH 8.1, 1 mM phenylmethylsulphonyl fluoride, 1 ^g/ml leupeptin, and 
1 jug/ml aprotinin). The resulting chromatin solution was sonicated to generate 
300-1 000 bp DNA fragments. After microcentrifugation, the supernatant was 
diluted with dilution buffer (0.01% sodium dodecyl sulphate, 1.1% Triton X-100, 
1.2 mM EDTA, 16.7mM Tris-chloride pH 8.1, 167 mM NaCI, 1 mM phenylmethyl- 
sulphonyl fluoride, 1 ^g/ml leupeptin, and 1 ^g/ml aprotinin) and precleared with 
protein A-Agarose/salmon sperm DNA and divided into aliquots. Five micrograms 
of p53 and HIF1a antibodies (Santa Cruz, Biotechnologies) or any antibody (as 
negative control) was added to the chromatin solution and incubated on a rotating 
platform overnight at 4°C. In parallel, we added acetylated H4-histone antibody 
(Millipore) or any antibody (as negative control) to a small fraction of the chromatin 
solution. Antibody-protein-DNA complexes were precipitated with protein 
A-Agarose/salmon sperm DNA. After centrifugation, the beads were washed 



and the protein-DNA complexes were eluted with 1% sodium dodecyl sulphate, 
100mM sodium carbonate. DNA-protein crosslinks were reversed by heating at 
65°C for 4h and the DNA was phenol extracted and ethanol precipitated. DNA 
fragments were analysed by real-time q-PCR using SybrGreen and primers 
specific for the CDKN1A (p21), as p53 activation-positive control, and PGC-1/?. 

p53 PCR amplification and sequencing. All 10 coding exons of TP53 
and exon-intron boundaries of the gene were amplified using Fast Start Taq DNA 
Polymerase (Roche Diagnostics). PCR amplification was performed in a final 
volume of 25 fi\ in a 9700 thermal cycler (Life Technologies), and PCR products 
were purified with Multiscreen PCR clean-up Filter Plates (Millipore). Primer 
sequences and PCR conditions for all reactions are available on request. 
Sequencing was performed with BigDye v1 .1 (Life Technologies) on both strands 
according to manufacturer's instructions. Briefly, 0.32 /uU of primer was used in 
10^1 final reaction, together with 0.5^1 of BigDyevI .1 and 2 [A of Buffer (2 x ). 
Upon cycling (4-min elongation) the amplicons were precipitated and run in an 
ABI3730 Genetic Analyzer automated sequencing machine (Life Technologies). 
Electropherograms were analysed using Chromas Lite v.2.01 (Technelysium, 
South Brisbane, QLD, Australia). 

TP53, MDM2, HIF1a, and HIF1a cloning, transfection, and 
degradation evaluation. TP53, MDM2, and HIF1a full-length cDNAs were 
cloned into pcDNA3.1 -vector (Life Technologies). To generate a form of HIF1a 
that is not subjected to hydroxylases degradation, HIF1a, and is constitutively 
active, amino acids P402, P564, and P803 of HIF1a were replaced with alanine 
residues using QuikChange Site-Directed Mutagenesis Kits (Agilent, Santa Clara, 
CA, USA). Cells were transfected using X-tremeGENE HP DNA Transfection 
Reagent (Roche Diagnostics), as described by the manufacturer. Mutagenesis 
primers are available on request. Treatments with chx (50 /^g/ml) were performed 
at 24 h after MDM2 transfection and cells collected after 0, 10, and 20 min of 
incubation. 

Senescence-associated /?-Gal assay. Increase in endogenous /?-Gal 
activity at higher pH than optimal are unequivocal markers of the switch towards a 
senescent phenotype. 55 Increased /?-Gal activity was determined as previously 
described. 55 The assay was performed on irradiated cells after recovery. Cells 
were washed twice with PBS and then incubated for 7 min with a fixation solution 
(2% formaldehyde and 0.2% glutaraldehyde in PBS). Cells were washed five times 
with PBS and incubated with a staining solution (40 mM citric acid/Na phosphate 
buffer pH 6, 1 mg/ml X-gal, 5mM K 4 [Fe(CN) 6 ]3H 2 0, 5mM K 3 [Fe(CN) 6 ], 150mM 
NaCI and 1 mM MgCI 2 in ddH 2 0) at 37°C until /?-Gal staining became visible. Cells 
were finally washed five times with PBS. 

Rectal biopsies collection. The study included 32 male patients, 18 years 
of age and older with prostate cancer. All patients underwent three-dimensional 
conformal radiotherapy and, by using 4-5-field techniques, a minimum energy of 
6MV was applied on linear accelerators. A computed tomography scan was 
performed before radiotherapy to enhance the reliability of the rectal dose-volume 
histogram. Dosing schedules followed the institutional protocols. During the 
treatment, patients laid in the supine position and were immobilized at the pelvis. 
Over 6- to 7-week period, using daily fractionation of 2Gy for 5 days a week, 
patients received a total dose of external-beam radiation ranging from 66 to 74 Gy, 
given in 33-37 fractions. The radiation dose was prescribed to the isocentre. 56 
During the examination by rectosigmoidoscopy, bioptic sampling was performed 
from the anterior wall of the rectum, 10 cm above the anal verge, before and 1 
month after the radiotherapy. Samples were snap frozen and kept at - 80°C until 
the molecular analysis was performed. Patients were enrolled following internal 
ethical committee procedures. 

Statistical analysis. Student's f-test was used for comparing two groups; 
ANOVA test was used for comparing more than two groups. Paired f-test was used 
for comparing the mtDNA copy number in pre- and post-IR single biopsies and 
Wilcoxon Rank tests were used for median comparison of all pre- and post-IR 
biopsies. P-value <0.05 was considered for any statistically significant association. 
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